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Abstract; Pseudodendrothrips mori is a serious pest to mulberry trees and its infestation affects the qualitative and 
quantitative characteristics of mulberry leaves consumed by silkworm. The spatial distribution of mulberry thrips within 
plant and orchard was assessed using Taylor’s power law and Morisita’s index of dispersion. The results showed that 
P . mori populations were localized in certain parts of the trees and fields. Distribution of thrips within plant showed that 
the larval instars were located on lower canopy (leaves 5- 10 from top ) of trees but adults preferred the upper canopy 
(leaves 1-5 from top ). The thrips density did not vary significantly among the leaf direction inside the trees. The thrips 
tended to be more accumulated on the trees located in East, South and North than in Center and West of orchard. The dis- 
tributions of P . mori adults and larvae on the leaves were aggregated . 
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The sericulture has a history of over 4000 years in low development and increases larval mortality; conse- 
east Asia and is an important part of the national econ- quently, there is a drop in the production of silkworm 
omy. Mulberry leaves are the exclusive food for the cocoons and a reduction in profit ( Etebari et al, 2004 ). 
silkworm, Bombyx mori. Feeding silkworms with mul- P . mori was recorded for the first time in 1999 in 


berry leaves that have been damaged by mulberry mulberry orchards of Guilan province, North of Iran 
thrips, Pseudodendrothrips mori ( Niwa, 1908 ) causes (Etebari et al, 1999). It was reported that feeding of 
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thrips decreases the leaf protein and moisture contents 
by 17.8 and 3.57%, respectively. Other investiga- 
tions demonstrated that the life cycle of this insect took 
29.95 days from egg to adult under laboratory condi- 
tions and had four generations in a year in Guilan 
province of Iran ( Jalali et al, 2001; Etebari et al, 
1999 ). All stages caused damages to trees by sap ex- 
traction, which resulted in leaf necrosis, drying of 
growing shoots and droping of premature leaves ( Cap- 
pellozza, 1987). Laboratory and field investigations 
showed that Kemuchi variety was the most susceptible 
variety and Shin-Ichinose was relatively resistant to this 
insect in Guilan ( Etebari et al, 2000 ). 

The importance of spatial heterogeneity comes 
from its main role in ecological theories and its practi- 
cal role in population sampling theory and in develop- 
ment of rational pests, management strategies as well 
as enhancement and conservation of natural enemies 
(Wang et al, 2004). The first step in the development 
of sequential sampling programs is to determine the 
spatial distribution of population being sampled. This 
determination is necessary to select the appropriate pro- 
cedure for the development of sequential sampling sys- 
tem (Hamilton et al, 1998). Insect populations are 
known to the follow three distinct distribution patterns 
(Southwood, 1978), viz. (i) binomial, uniform or 
regular, ( ii ) random or poisson, and ( iii ) negative bi- 
nomial, aggregated or clumped. However there is no 
published information on within-plant and within-field 
distribution of P.mori. If thrips are highly aggregat- 
ed, more samples will be needed than their regular dis- 
tribution. The degree of aggregation may also reveal in- 
formation about thrips behaviour (Kirk, 1997). So, 
knowledge of the spatial distribution of P.mori is 
needed to provide information for the possible manage- 
ment strategy. Also economic threshold needs a careful 
sampling consideration for both thrips distribution with- 
in the tree and the among different trees in mulberry 
orchards. The present research hoped to point out the 
spatial distribution of this thrips in orchard that can be 
used in developing a sampling plan with an acceptable 


precise level. 
1 Materials and Methods 


1.1 Study site and Sampling 

The study was carried out in mulberry orchards of 
Iran Sericultural Research Center ( ISRC ), in Pasikhan 
region (North of Iran, Guilan Province ). In the study 


site, five individual trees were selected in each orienta- 
tion (North, South, West and East ) and also in the 
center of orchards. The sampling is started in July 20 
and leaves were collected every week until September 
20, also the pervious studies showed that at this time, 
the population of thrips is in its highest density ( Jalai 
et al, 2001 ). Each tree was divided into upper (leaves 
1-5 from top ) and lower canopy (leaves 5- 10 from 
top ) and in each canopy four leaves from different ori- 
entations (N, S, W, and E ) were collected separately 
in plastic bags and kept in fridge for 2 hours. The 
number of first, second instar larvae and adults of 
P . mori was counted on leaves. 

1.2 Statistical and Data analysis 

The spatial distribution patterns of P . mori were 
determined by using Taylor’s power law ( Taylor, 
1961). Taylor’s power law parameters were obtained 
by the regression of logyo-transformed variances, s”, on 
logio-transformed mean numbers of P . mori adults and 
larvae per sample, i.e., by means of the linear regres- 
sion model: log s* = log a + b log x (Taylor 1961). 
According to this model b value > 1 denotes a popula- 
tion with an aggregated distribution, 6 value signifi- 
cantly <1 denotes a regular distribution, and 6 value 
not significantly different from 1 denotes a random dis- 
tribution. The fit of each data set to the linear regres- 
sion model was evaluated by calculating the r? value. 
The Student t- test was used to determine if the slopes 
(b values ) obtained by means of the linear regression 
procedure were equal to 1, significantly greater than 1 
or significantly less than 1. 

Also, Index of Dispersion (ID ) was calculated for 
all series of data with the formulae ID = S?/m, where 
if ID = 1 the dispersion is random, ID = 0 shows a uni- 
form distribution and for clumped pattern ID should be 
more than 1. In order to compare different models MI 
(Morisita’s Index of Dispersion ) and MSI ( Morisita’s 
Standardized Index of Dispersion ) indices were mea- 
sured to check if MSI is more than zero the clumped 
situation or less than zero the uniform situation is 
shown ( Pourbabaei, 2004 ). 

Weekly thrips counts were analyzed to determine 
significant level of thrips density in the sampling parts 
ANOVA(SAS, 1997). For the analysis, data were 
transformed into natural logarithms, but all the data 
are presented in original units as means + SE. Means 
were separated using Duncan Multiple Range Test 


(DMRT ). 
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2 Results and Discussion 


The obtained results showed significant differences 
(P<0.05, df=1) among mulberry thrips population 
collected at two different heights (lower and upper 
canopy ), five positions of trees in orchard ( North, 
South, East, West and center ), three different stages 
of thrips development (1st and 2nd larva and adult ) 
and three varieties of mulberry. But there was no sig- 
nificant difference among thrips populations at four di- 
rections of the trees (Tab. 1). 


Tab. 1 Analysis of variances for characters and interactions 
for mulberry thrips, Pseudodendrothrips mori 
Parameters df F P 
Position of tree 4 6.01 0.0001 * 
iL 95.96 < 0.0001* 
2 64.39 <0.0001* 
3 1.03 0.38 
Variety 2 64.39 0.0001 * 
Height of canopy x Leaf direction 3 0.86 0.46 
2 
2. 
6 
4 








Height of canopy 
Thrips instar 
Leaf direction 


Height of canopy x Thrips instar 58.53 0.0001* 
Height of canopy x Variety 40.46 0.0001* 
Thrips instar x Leaf direction 0.70 0.65 

Thrips instar x Variety 31.27 0.0001 * 


Variety x Leaf direction 6 1.22 0.29 
“If P <0.05, then the means are significantly different 


Abundance of mulberry thrips was higher in East, 
South and North trees in orchard than Center and West 
trees especially in Km variety (Fig. 1). Pearsall and 
Myers ( 2001 ) showed that flight of Western flower 
thrips, Frankliniella occidentalis ( Pergnade ) in Nec- 
tarine orchards of British Columbia, appeared likely to 
be wind driven into some orchards. Also it was report- 
ed that western flower thrips densities were higher on 
sticky cards in orchards located in the South comparing 
with the North ( Pearsall and Myers, 2001 ). From this 
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Fig. 1 Mulberry thrips density in five positions of 


different varieties in orchard 


research it is suggested that the average thrips densities 
in all varieties were higher on Eastern (Mean + SE: 
7.553 + 18.11), Southern ( 7.091 + 24.11 ) and 
Northern ( 5.954 + 20.966 ) parts of trees. Although 
the thrips population on Eastern and Southern parts of 
trees were higher than North, these differences were 
statistically not significant ( P > 0.05, df=4). 

In all mulberry varieties, larval instars had the 
highest density on the lower canopy, while adults pre- 
ferred the upper parts of the plant canopy ( Fig. 2). 
Kirk (1997 ) reported that the site occupied by thrips 
can differ between adult and larvae, in which adults are 
often nearer the top of the plant and larvae lower down. 
At finer scale, there can be differences between leaf 
ages and between the upper and lower sides of leaves. It 
suggested that the distribution patterns are mainly relat- 
ed to feeding behavior but oviposition and mating behav- 
ior have an important role in thrips distribution . 

Pearsall and Myers ( 2001 ) showed that western 
flower thrips moved into orchards equally at ground lev- 
el and at the 1 m height, which was the height of the 
tree canopy. This might be due to the tendency of 
avoiding the effects of wind at higher levels of ground 
and tree for the thrips. Temperature is another factor 
known to influence aerial dispersal of invertebrates 
(Taylor, 1963 ). That could be responsible for estab- 
lishment of thrips in different tree height levels. Since 
adult thrips are thought to make flights just above the 
plant canopy ( Brodsgaard, 1989 ), it is reasonable to 
expect that most adults would be found in the upper 
canopy. The limited distribution of immatures within 
the lower canopy suggests that P.mori may prefer a 
protected habitat from direct sunlight and wind, as do 
other species of thrips ( Picket et al, 1988). In this 
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Fig. 2 Mulberry thrips density in three different develop- 
mental stages at two heights ( lower and upper 


canopy ) in three varieties 
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case, Seal et al (2006) found that Scirtothrips dorsalis 
Hood adults and larvae were abundant on terminal 
leaves of “Scotch Bonnet’ pepper plants followed by 
middle leaves, lower leaves, and reproductive parts, 
while Reitz ( 2002 ) mentioned that adults of Franklin- 
iella species tended to be in flowers in the upper part of 
tomatoes whereas more immature thrips occur in the 
lower parts of the plant canopy. 

The analysis of thrips densities on the leaves in 
the four directions of the trees showed that there was no 
difference between directions of leaves in the trees 
(Fig. 3)( P =0.38, df =3). Also interactions of leaf 
direction with other variations were not significantly dif- 


ferent (Tab. 1). 
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Fig. 3 Mulberry thrips density in four directions of 
each variety of mulberry trees 


The results showed that the thrips was available in 
larger numbers in KM variety (Mean + SE: 11.583 + 
28.434) than the two other varieties (I: 2.5685 + 
8.784, SI: 1.9789 + 7.319). Etebari et al (2000 ) 
with field and laboratory experiments revealed that KM 
was the most susceptible and SI was the most resistant 
to P . mori. Underwood et al ( 2005 ) suggested that in- 
duced and constitutive resistance could affect the spa- 
tial distribution of herbivore damage differently. 

In recent study, the most abundant thrips stage 
was the Ist instar larva (mean: 12.0617 ) in compari- 
son of the 2nd instar (mean: 2.2724) and adult 
(mean: 1.7968), as has been noted by other re- 
searchers ( Jalali et al, 2001 ). It may because of pres- 
ence and efficiency of the thrips’ natural enemies on 
2nd instar larva. 

In all of the thrips life stages ( Ist and 2nd larva 
and adult ) b value was higher than Ist, hence the 
thrips were aggregated on the trees (Tab. 2). The 
highest slope value was noted for the adult ( b = 


6.4281). Also the values of r? obtained with Taylor’s 
power law were higher for all stages. This indicates a 
good fit of model to the data on thrips using leaves as 
the sampling unit. Because thrips tend to aggregate in 
certain location for mating ( Reitz, 2002), maybe it is 
reasonable to explain significant difference of adults ag- 
gregation in comparison of larval instars in Tab. 2. 
This aggregated pattern of distribution in field is also 
typical for other thrips species such as Scirtothrips dor- 
salis (Seal et al, 2006), Thrips palmi (Seal 1996, 
Seal and Stansly, 2000) and Frankliniella occidentalis 
(Rhainds and Shipp, 2003 ). 


Tab. 2 Distribution of mulberry thrips based on 
Taylor ’s Power Law 








Stage re a b 

Adult 0.803 6.473 6.428 AGG* 
Ist larva 0.954 0.182 1.947 AGG 
2nd larva 0.944 1:35: 2.774 AGG 





* AGG: b> 1 shows aggregated distribution. 


The data in Tab. 3, which outline MSI, MI, ID 
indices also shows the aggregate dispersion in all the 
larval stages and adults. ID is not unique in all cases. 
The first instar larvae has more indics and this demon- 
strates aggregate dispersion in this stage of insect life 
(Tab. 3). Therefore, all models presented an aggre- 
gate dispersion for mulberry thrips. 

Understanding the causes and consequences of 
dispersal, it is crucial to determine the impact of dis- 
persal on the spatial structure of animal populations 
(Nathan 2001 ). In several phytophageous insects, 
crowding of the host plant is correlated with a high in- 
cidence of dispersal, either as a result of high popula- 
tion density per se or intense feeding by herbivores on 
their host plant (Denno and Peterson 1995, Herzig 
1995, Herzig and Root 1996, Dixon 1998, Rhainds et 
al, 2002}. Density-dependent dispersal has a strong 
influence on the population dynamics of phytophageous 
insects because it stabilizes local populations and rela- 
tively generates uniform, synchronous distribution of 
insects, especially when adults preferentially colonize 
hosts with a low abundance of conspecies ( Ruxton and 
Rohani 1998 ). 

Spatial distribution of insects within trees may be 
affected by numerous factors including parasitism and 
predatism ( Hassell, 1978; Ehler, 1986), including 
dispersal ability of immatures or adults ( Witham, 
1981), competition ( Hails and Crawley, 1991), host 
plant resistance ( Underwood et al, 2005) and plant 
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Tab. 3 Some statistical parameters and distribution index for mulberry thrips 











Lower Crown Upper Crown 
TE Statistical F r k md 
Dregont., Shanes Aduk : 7 Adult ! ? 
Instar Instar Instar Instar 
Northern Mean 1.600 25.22 4.367 1.978 2.522 0.222 
SE 0.218 4.630 1.054 0.332 0.364 0.070 
ID 2.680 76.49 22.87 5.010 4.733 2.000 
MI 2.045 3.961 5.967 3.017 2.470 5.684 
MSI 0.505 0.516 0.527 0.510 0.508 0.518 
Southern Mean 1.778 20.85 3.456 2.089 2.933 0.478 
SE 0.310 3.900 0.751 0.350 0.772 0.220 
ID 4.876 65.65 13.92 5.268 18.29 9.135 
MI 3.170 4.070 4.711 3.031 6.852 18.23 
MSI 0.511 0.516 0.520 0.511 0.532 0.593 
Eastern Mean 1.633 16.28 3.376 1.956 2.089 0.289 
SE 0.289 2.889 0.636 0.366 0.392 0.070 
1D 4.615 46.12 10.823 6.181 6.634 1.886 
MI 3.204 3.741 3.895 3.635 3.681 4.154 
MSI 0.511 0.515 0.516 0.514 0.514 0.512 
Western Mean 1.533 23.06 5.689 2.178 2.311 0.289 
SE 0.242 4.607 1.864 0.309 0.555 0.111 
1D 3.447 82.82 54.97 3.958 11.98 3.831 
MI 2.590 4.509 10.401 2.350 5.723 11.07 
MSI 0.508 0.520 0.553 0.507 0.525 0.551 


SE: Standard Error; ID: Index of Distribution; MI; Morista’s Index of Distribution; MSI: Morista Standardized Index of Distribution 


phenology ( Pearsall and Myers, 2001; Faeth et al, 
1981). These results indicate that the mulberry thrips 
is more abundant on KM variety and in Ist instar lar- 
va. P.mori is found mainly on the lower leaves of 
mulberry, which may avoid the influences of wind 
speeds and high temperature on the thrips. Also the 
thrips was more abundant on trees located in East, 
South and Northern parts of orchard. It suggests that 
the sampling should involve all the surface of orchards. 
The thrips distribution was not related to leaf direction 
in the tree. 

This study provides the first experiment on the 
spatial distribution of P . mori. Our data indicate that 
mulberry thrips were spatially distributed on the trees. 


Also we found that P . mori populations were localized 
References: 


Brodsgaard HF. 1989. Colored sticky traps for Frankliniella occidentalis 
(Thysanoptera: Thripidae ) in glasshouses [ J ]. Journal of Applied 
Entomology , 107: 136-140. 

Cappellozza L, Miotto F. 1987. Observations on the bio-ethological be- 
haviour of Pseudodendrothrips mori ( Niwa )( Thys.: Terebrantia ) in 
some zones of Italy [J]. Redia , 70: 523-533. 

Denno RF, Peterson MA. 1995. Density-dependent dispersal and its 
consequences for population dynamics [ A ]. In: Cappuccino N, 
Price PW (eds. ), Population Dynamics: New Approaches and Syn- 
thesis [M]. San Diego, CA: Academic press. 113-130. 

Dixon AFG. 1998. Aphid Ecology: An Optimization Approach [ M ]. 
London, United Kingdom: Chapman & Hall. 


in certain parts of the tree. The thrips is often recorded 
with high population densities in mulberry orchards in 
Iran and elsewhere ( Cappellozza and Miotto, 1987; 
Etebari et al, 1999; Etebari atal, 2004; Jalali et al, 
2001; Yin et al, 1994). Results in this paper can be 
used for the development of sampling protocols and also 
in the management of P . mori in orchards. Additional 
research is required in order to assess these variables 
for a wider range of cases, probably the effect of para- 
sitism and predation on this thrips and finally to desig- 
nate an action threshold for this pest in this region. 


Acknowledgments: This work was financed by 
the University of Guilan and we thank L. Matindoost 
and A. R. Bizhannia for their technical assistance . 


Ehler LE. 1986. Distribution of progeny in two ineffective parasites of a 
gall midge ( Dipt.: Cecidomyiidae )[ J ]. Environmental Entomolo- 
gy, 15: 1268-1271. 

Etebari K, Jalali J, Taksokhan M. 1999. First record of mulberry thrips 
Pseudodendrothrips mori NIWA( Thy . : Thripidae ) from mulberry or- 
chards of north Iran [J]. Journal of Entomological Society of Iran , 
18 (1/2): 45. 

Etebari K, Jalali J, Taksokhan M. 2000. Investigation on host preference 
of Pseudodendrothrips mori Niwa (Thy; Thripidae ) for varieties of 
mulberry in the field and laboratory [J]. Iranian Journal of Agricul- 
tural Science , 31 (3 ): 489-495. 

Etebari K, Matindoost L, Singh RN. 2004. Decision tools for mulberry 


270 Zoological Research 


Vol. 28 





thrips Pseudodendrathrips mori ( Niwa, 1908 ) management in seri- 
cultural regions: An overvie[ J]. Entomologia Sinica , 11 (4 ): 243 
=255. 

Faeth ST, Mopper S, Simberloff D. 1981. Abundance and diversity of 
leaf-minig insects on three oak host species: Effects of host plant 
phenology and nitrogen content of leaves [J]. Oikos, 37: 238-251. 

Hails RS, Crawley MJ. 1991. The population dynamics of an alien in- 
sect: Andricus quercuscalicis (Hym.: Cynipidae )[ J]. Journal of 
Animal Ecology , 60: 545-562. 

Hassell MP. 1985. Insect natural enemies as regulating factors [ J J. 
Journal of Animal Ecology, 54: 323-334. 

Herzigal 1995. Effects of population density on longdistance dispersal in 
the goldenrod beetle [J]. Trirhabda virgata. Ecology, 76: 2044- 
2054. 

Herzig AL, Root RB. 1996. Colonization of host patches following long- 
distance dispersal by a goldenrod beetle [ J ]. Trirhabda virgata. E- 
cological Entomology, 21: 344-351. 

Jalali sendi J, Hasheminia M, Etebari K, Shojaee M. 2001. The labora- 
tory biology and seasonal occurrence of mulberry thrips Pseudoden- 
drothrips mori Niwa ( Thys. , Thripidae ) in Guilan Province, ( Pas- 
sikhan region )[ J]. Journal of Entomological Society of lran, 21 
(2): 87-100. 

Kakimoto K, Inoue H, Yoshida T. 2002. Seasonal changes in population 
density and species composition of Orius spp. in mulberry fields 
[J]. Japanese Journal of Applied Entomology and Zoology , 46 ( 4 ): 
209-215. 

Kirk WDJ. 1997. Distribution, Abundance and Population Dynamics 
[A]. In: Lewis T (eds. ), Thrips as Crop Pest [M ]. CAB Interna- 
tional, 217-258. 

Nathan R. 2001. The challenges of studying dispersal [ J]. Trends in E- 
cology and Evolution , 16: 481-483. 

Pickett CH, Wilson LT, Gonzalez D. 1988. Population dynamics and 
within-plant distribution of the western flower thrips ( Thysanoptera : 
Thripidae ), an early-season predator of spider mites infesting cotton 
[J]. Environmental Entomology, 17 (3): 551-559. 

Pearsall IA, Myers JH. 2001. Spatial and temporal patterns of dispersal 
of Western flower thrips ( Thysanoptera: Thripidae ) in nectarine or- 
chards in British Columbia [J]. Journal of Economic Entomology , 
94 (4): 831-843. 

Pourbabaei H. 2004. Statistical Ecology [ M ]. The University of Guilan. 
428 p. 

Prasad KS, Sreedhar $, Singhvi NR, Kodandaramaiah J, Sen AK. 
2002. Post thrips infestation biochemical changes in leaves of mul- 


berry ( Morus spp. )[J]. Plant Archives, 2 (1): 85-88. 

Reitz SR. 2002. Seasonal and within plant distribution of Frankliniella 
thrips ( Thysanoptera: Thripidae ) in North Florida tomatoes [ J ]. 
Florida Entomologist , 85 (3 ): 431-439. 

Rhainds M, Gries G, Ho CT, Chew PS. 2002. Dispersal by bagworm 
larvae, Metisa plana: Effects ofpopulation density, larval sex, and 
host plant attributes [J]. Ecological Entomology , 27: 204-212. 

Rhainds M, Shipp L. 2003. Dispersal of adult Western Flower Thrips 
(Thysanoptera: Thripidae ) on Chrysanthemum Plants: Impact of 
feeding-induced senescence of inflorescences [ J ]. Environmental 
Entomology, 32 (5): 1056-1065. 

Ruxton GD, Rohani P. 1998. Fitmess-dependent dispersal in meta popu- 
lations and its consequences for persistence and synchrony [ J ]. 
Journal of Applied Ecology, 67: 530-539. 

SAS Institute. 1995. SAS/STAT user’s guide, release 6.0[ CP ]. Cary, 
NC: SAS Institute. 

Seal DR, Ciomperlik M, Richard ML, Klassen W. 2006. Distribution of 
the Chilli thrips Scirtothrips dorsalis ( Thysanoptera: Thripidae ) 
within pepper plants and within pepper fields on St. Vincent [J ]. 
Florida Entomologist , 89 (3): 311-320. 

Song ZY. 1981. A discussion on the method of monitoring and forecasting 
mulberry thrips [J]. Jiangsu Sericulture, 2: 27. (in Chinese ) 
Southwood TRE. 1978. Ecological Methods with Particular Reference to 

the Study of Insect Populations [M ]. New York: Wiley. 

Taylor LR. 1961. Aggregation, variance and the mean [J]. Nature 
( London ), 189: 732—735. 

Taylor LR. 1984. Assessing and interpreting the spatial distributions of in- 
sect populations [J ]. Annual Review of Entomology, 29: 321-357. 

Taylor LR. 1963. Analysis of the effect of temperature on insects in flight 
[J]. Journal of Animal Ecology , 32: 99-117. 

Underwood N, Anderson K, Inouye BD. 2005. Induced vs. constitutive 
resistance and the spatial distribution of insect herbivore among 
plants [J]. Ecology, 86 (3): 594-602. 

Wang Z, Yuan Z, Song D, Zhu M. 2004. Analysis on the dynamics of 
spatial distribution pattern of mixed spider population in rice field 
[J]. Entomologia Sinica , 11 (4): 257-265. 

Whitham TG. 1981. Individual trees as heterogeneous environments: 
Adaptation to herbivory or epigenetic noise? [ A ]. In: Denno RF, 
Dingle H( eds. ), Insect Life Listory Pattern: Habitat and Geograph- 
ic Variation [M ]. New York: Academic Press. 9-27. 

Xi TZ, Zhu J. 1991. The regular pattern of occurrence of mulberry thrips 
and the control strategies [J]. Sichuan Sericulture, 2: 24-25. (in 
Chinese ) 


